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ABSTRACT 


A  computer  program  has  been  developed  as  a  result  of  this  project.  For  input 
data  consisting  of  the  spectral  distribution  of  the  incident  energy,  the  composition  of  the 
scatterer,  a, ’id  the  size  distribution  of  the  scatterer,  this  program  computes  the 
volumetric  scattering  intensity.  Preliminary  results  obtained  with  the  program  indi¬ 
cate  a  very  prominent  degree  of  lateral  scattering. 


INTRODUCTION 


To  afford  the  reader  a  fuller  understanding  of  this  project,  the  application  of 
radiant  energy  to  fire  spread,  which  is  actually  the  main  stimulus  for  this  study,  will 
initially  be  described.  Following  this  description  of  fire  spread,  the  history  of  the 
development  of  this  project  will  be  briefly  chronicled,  so  that  the  reader  may  more 
fully  appreciate  why  one  particular  small  segment  of  the  problem  was  taken  out  for 
study.  Finally,  somewhat  in  the  nature  of  a  summary,  the  history  of  the  Mie  problem 
will  be  considered. 


SUMMARY 


The  accomplishments  of  this  project  have  been: 

(1)  The  preparation  of  a  computer  code  by  which  numerical  data  for  the 
volumetric  scattering  intensity  can  be  obtained.  Incorpoi  ated  within  * 
the  program  is  a  plotting  option  by  which  the  data  may  be  portrayed 
graphically. 

(2)  Limited  experience  with  the  computer  program  has  indicated  tentative 
bounds  for  the  convergence  of  the  infinite  summation,  and  has  identi¬ 
fied  the  computational  steps  that  must  be  checked  to  ensure  convergence. 

(3)  The  output  data  obtained  to  date,  presented  as  a  plot  of  intensity  against 
scattering  angle,  demonstrate  the  very  strong  degree  of  lateral 
scattering  inherent  in  the  Mie  problem.  Successful  enhancement  of 
fire  spread  will  require  serious  consideration  to  both  the  gross  cov 
figuration  of  the  reflective  cloud,  and  the  minute  aspects  of  the 
scatterer  —  composition  and  size. 
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CONCLUSIONS 


As  a  result  of  this  program,  the  following  conclusions  were  reached: 

(1)  The  Mie  problem  can  be  programmed  with  nominal  difficulty,  and  for 
a  computer  operating  interval  of  three  minutes  maximum  will  produce 
a  curve  of  scattered  intensity  against  scattering  angle  for  a  specified 
wavelength  and  scattering  particle.  If  the  problem  is  set  up  for  multiple 
inputs,  the  running  time  ir,  lengthened  about  one -half  minute  j  -  curve. 

(2)  The  versatility  of  the  program  should  be  capabl?  of  considerable  ex¬ 
tension;  however,  this  is  more  of  an  engineering  problem  than  a 

p:  og ramming  problem.  A  clearer  definition  of  the  fire-spread 
problem  would  stimulate  improvements  in  the  program. 


RECOMMENDATIONS 


Inherently  an  exploratory  program  such  as  this  would  allow  the  compilation  of 
a  rather  large  list  of  recommendations,  and  mos*  of  these  recommendations  would  in¬ 
volve  additional  work  to  extend  a  particular  aspect  of  the  study.  However,  3ince  there 
is  no  indication  that  an  extension  of  this  program  is  likely,  only  a  few  very  general 
recommendations  will  be  offered.  For  instance,  optical  data  for  the  various  candidate 
materials  to  be  used  for  reflective  clouds  are  not  available  in  a  form  readily  usable  by 
the  cloud  designer.  Some  efforts  leading  to  a  more  handy  presentation  would  seem 
justified.  Regarding  the  efforts  here,  several  aspects  of  the  computer  program  could 
be  modified  to  extend  the  range  and  versatility  of  this  program. 


BACKGROUND 


Application  of  Radiation  to  Fire  Spread 


In  any  analysis  of  mass  fires  as  either  an  offensive  or  a  defensive  weapon, 
serious  consideration  must  be  given  to  ignition,  the  initiation  of  the  fire,  and  to  propa¬ 
gation  of  the  fire,  at  least  during  the  early  stages.  This  cor  dideration  of  ignition  and 
propagation  mav  be  combined  into  the  general  term  "fire  spread".  The  use  of  mass 
fire  as  a  weapon  implies  enhanced  fire  spread.  A  discussion  of  fire  spread  necessarily 
involves  the  energy  balance  associated  with  the  fire.  One  very  significant  term  in  this 
energy  balance  is  the  radiant  energy  term,  for  radiant  energy  fills  the  role  of  a 
precursor  in  mass  fire  spread. 

An  early  realization  in  any  consideration  of  this  radiant  energy  term  is  the  fact 
that  much,  if  not  most,  of  the  radiant  energy  associated  with  an  open  fire  escapes 
skyward  and  is  lost  to  the  energy  balance.  It  is  then  rather  obvious  that  in  any  effort 
to  enhance  fire  spread,  an  attempt  should  be  made  to  return  earthward  a  significant 
portion  of  this  skyward-bound  radiant  energy.  In  a  military  situation,  the  attempted 
ignition  will  utilize  incandiary  bombs  and  here  too  much,  if  not  most,  of  the  radiant 
energy  is  lost  skyward. 
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In  elementary  terms,  the  escaping  radiant  energy  may  be  returned  earthward  by 
spreading  a  reflective  blanket  above  the  fire  zone.  In  military  terms,  this  reflective 
blanket  can  be  produced  by  spreading  a  reflective  smoke  screen  or  a  reflective  fog 
above  the  fire.  Snv'ke  generators,  or  fog  generators,  either  earthbound  or  airborne, 
represent  a  reasonably  well-developed  technology,  so  suitable  delivery  means  are 
within  the  current  state  of  the  art.  The  only  remaining  question  is  the  effectiveness  of 
the  smoke  cloud  in  returning  energy  earthward.  It  is  to  this  question  that  the  efforts 
represented  by  this  report  were  directed. 


History  of  RACIC/ LWL  Interchange 


In  September  of  1965,  a  request  was  made  by  the  Limited  War  Laboratory  at 
Aberdeen  Proving  Ground  to  Battelle's  RACIC- Washington  office  for  information  related 
to  the  reflective  properties  of  smoke  clouds.  The  response  to  this  request  initially  took 
the  form  of  a  search  of  the  RACIC- Columbus  files.  To  supplement  the  RACIC  holdings, 
a  number  of  inqux.  ies  were  made  to  investigators  whose  general  scope  of  study  sug¬ 
gested  that  they  might  be  able  to  contribute  pertinent  information.  The  activities  re¬ 
lated  to  this  initial  request  were  concluded  in  late  October  1965  and  reported  in  a 
letter  to  LWL  tabulating  the  bibliographical  data  and  other  information  obtained. 

From  the  onset,  this  problem  was  of  a  strong  heuristic  nature  and  generated 
many  discussions  regarding  the  probable  implications.  Further,  the  amount  of  data 
gathered  during  the  literature  search  was  rather  small  and  of  quite  disappointing 
quantitative  content.  The  extended  needs  of  LWL  were  discussed  and  the  availability 
of  RACIC  assistance  was  made  known.  Subsequently,  LWL  requested  such  assistance 
from  RACIC,  and  in  late  January  of  1966  this  study  was  authorized. 

During  the  time  interval  in  which  authorization  of  the  project  was  considered, 
there  was  personal  communication  between  LWL  personnel  and  RACIC  personnel. 

During  these  conversations  the  broad  field  of  radiant  energy  in  fire  spread  was  sifted, 
to  give  one  aspect  which  might  be  attacked  with  the  resources  at  hand.  The  fin  1 
mutual  decision  was  to  prepare  a  computer  program  which  would  provide  an  estimate 
of  the  volumetric  scattering  efficiency  for  any  specified  cloud.  This  decision  was  sub¬ 
sequently  confirmed  by  correspondence*. 

Statement  of  the  LWL  Problem 


The  field  of  radiant  energy  and  its  interaction  with  matter  is  exceedingly  broad. 
Any  given  effort  can  at  best  attack  only  a  limited  region,  as  these  efforts  were  designed 
to  do.  To  enhance  appreciation  of  the  results  stated  in  this  report,  the  intended  purpose 
will  be  clearly  and  rather  restrictively  drawn.  One  element  in  tne  planning  of  future 
work  to  develop  a  newly  hypothesized  smoke  is  the  degree  of  expectation  that  thi  smoke 
would  achieve  the  results  postulated.  The  need  which  led  to  this  program  was  for  an 
established  methodology  which  would  quickly  give  a  working  estimate  of  the  reflective 
effectiveness  and  thus  assist  in  planning  a  development  program.  The  intent  here  was 
to  establish  the  methodology  to  give  the  working  estimate  by  programming  the  Mie 
solution  in  a  quite  restrictive  fashion. 


Reference  teletype  W-21  dated  December  17,  1965,  from  F,  J,  Barone /RACIC -W  to  J.  T.  Brown/R.  E.  Hess/RACIC-C,  and  letter 
CRD-AM-5A  dated  December  13,  1965,  from  Lt.  Col.  J.  F.  Mattes  on,  Chief,  Operations  and  Analysis  Division,  Limited  War 
Laboratory,  to  F.  j.  Barone/RACIC-W. 
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Description  of  the  Physical  Problem 


The  Mie  scattering  problem  is  the  description  of  the  interaction  of  a  plane 
electromagnetic  wave  with  a  spherical  scatterer.  This  problem  is  applicable  to  any 
portion  of  the  electromagnetic  spectrum  and  may  be  used  to  describe  the  scattering  of 
radar  waves,  of  thermal  radiation,  or  of  visible  light.  The  application  here  is  to 
thermal  radiation  as  it  might  be  scattered  by  a  cnemical  smoke  used  in  military 
operations. 

The  elementary  definition  of  electromagnetic  radiation  employs  an  electric 
vector  and  a  magnetic  vector  which  oscillate  perpendicularly  to  each  other  and  which 
define  a  pL  ae  whose  normal  is  the  direction  of  propagation  of  the  electromagnetic  wave. 
Strictly  this  description  fits  only  one  wave  train,  and  a  beam  of  radiation  would  consist 
of  many  such  wave  trains.  For  an  unpolarized  beam,  there  is  no  preferred  orientation 
so  that  the  electric  vectors,  hence  the  magnetic  field  also,  art  randomly  oriented 
with  all  directions  equally  probable.  Under  these  circumstances,  the  electric  field  and 
magnetic  vectors  can  be  each  represented  by  two  perpendicular  vectors  of  equal 
magnitude. 

If  this  electromagnetic  wave  traverses  an  electron,  the  electron  responds  to  the 
electric  field  of  the  wave  and  is  force  J  into  oscillation  by  <diis  field.  This  moving 
electron  now  gives  rise  to  a  magnetic  field  ant  since  the  movement  of  the  electron  is 
accelerated  motion,  the  resulting  magnetic  field  has  a  timewise  variation.  Since  a-> 
changing  magnetic  field  gives  rise  to  a  changing  electric  field,  the  motion  of  the  electron 
results  in  the  propagation  of  a  secondary  electromagnetic  wave.  The  electron,  in 
responding  to  the  incident  electric  field,  oscillates  in  the  plane  of  this  electric  field; 
the  field  generated  by  the  movement  of  the  electron  is  in  a  plane  perpendicular  to  the 
plane  of  oscillation.  Although  this  description  is  greatly  oversimplified,  it  serves  to 
demonstrate  that  the  interaction  of  an  incid.nt  field  and  a  responsive  electron  removes 
energy  from  the  original  beam  and  scatters  it  in  a  plane  perpendicular  to  the  initial 
beam. 


The  basic  task  now  is  to  relate  this  scattering  to  a  three-dimensional  geometry. 
For  the  very  brief  explanation  permitted  here,  the  earth  will  be  taken  as  a  flat  surface 
defining  the  xy-plane  and  the  perpendicular  to  this  plane  is  the  local  vertical  defined  as 
the  z-axis.  The  plane  electromagnetic  wave  under  discussion  will,  of  course,  be  un¬ 
polarized  and  will  be  described  by  vectors  Eox  and  E0y  of  equal  magnitude.  Associated 
v  \th  the  electric  vector  Eox  is  a  magnetic  vector  H0y  and  associated  with  the  electric 
vector  E0y  is  a  magnetic  vector  -Hox.  The  one  negative  vector  is  necessary  to  give 
the  same  propagation  direction  to  each  pair  of  component  vectors. 

When  this  plane  wave  passes  over  an  electron,  the  Eox  vector  causes  the  electron 
to  oscillate  in  the  xy-plane  along  the  x-direction  while  the  E0y  vector  causes  the 
electron  to  also  oscillate  in  the  xy-plane  but  along  the  y-direction.  The  electric  field 
variations  are  expressed  in  sinusoidal  terms,  so  the  motion  of  the  election  is  simple 
harmonic,  at  least  in  the  ideal  sense.  In  simple  harmonic  motion  the  moving  particle 
has  zero  velocity  at  the  ends  of  its  travel  and  a  maximum  veloc'ty  at  its  midpoint. 

The  strength  of  the  magnetic  field  created  by  a  current  is  proportional  to  the  current, 
and  the  flux  lines  are  circles  centered  on  the  conductor.  The  moving  electron  is  a 
current,  and  since  the  velocity  of  the  electron  is  changing  sinusoidally,  the  current  it 
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represents  is  also  changing  sinusoidally.  Thus,  the  field  strength  is  greatest  when 
the  electron  is  passing  through  its  midpoint  and  this  field  strength  decays  to  zero  at  the 
ends  of  its  travel. 

The  oscillations  of  the  electron  due  to  the  incident  electr;.c  vector  Eox  give  rise 
to  the  scattered  magnetic  vectors  Hsy  and  HSz.  By  Lenz's  Law  the  changing  magnetic 
field  resists  the  current  variation.  Thus,  there  is  a  scattered  electric  vector  -Esx 
associated  with  each  of  the  scattered  magnetic  vectors.  The  combination  of  vectors 
-Esx  and  HSy  has  as  its  direction  of  propagation  the  minus  z-direction;  in  other  words, 
this  combination  represents  the  radiation  scattered  along  —  both  forward  and  backward  — 
the  direction  of  the  incident  ray.  The  other  combination  of  vectors  -Esx  and  Hez  has 
as  its  direction  of  propagation  the  y-direction;  in  other  words,  this  combination  repre¬ 
sents  the  radiation  scattered  laterally. 

Similarly  the  oscillation  of  the  electron  due  to  the  vector  E0y  gives  rise  to 
scattered  magnetic  vectors  Hsx  and  Hsz,  and  associated  with  each  of  these  magnetic 
vectors  is  an  electric  vector  -ESy.  The  combination  of  vectors  -ESy  and  Hsx  has  as 
its  direction  of  propagation  the  z-direction,  hence  it  represents  the  radiation  scattered 
coincidentally  with  the  incident  ray.  The  combination  of  vectors  -E8V  and  Hsz  has  as 
its  direction  of  propagation  the  minus  x-direction,  so  it  represents  the  radiation 
scattered  laterally. 

Since  the  electron  is  in  simple  harmonic  motion  along  the  two  axes  in  its  plane  of 
oscillation,  there  is  symmetry  about  the  midpoint  of  its  excursion.  Thus,  the  energy 
scattered  along  the  direction  of  the  incident  ray  is  scattered  equally  in  the  earthward 
and  skyward  directions.  Similarly  for  the  laterally  scattered  energy,  equal  amounts 
are  scattered  to  the  right  and  to  the  left  and  equal  amounts  forward  and  aft.  The 
accounting  of  this  scattered  energy  is  the  task  of  this  propram. 


History  of  the  Mie  Solution 


Once  the  mathematical  formalism  related  to  the  electromagnetic  field  had  been 
developed  by  Maxwell  in  the  middle  of  the  nineteenth  century,  many  attempts  were 
made  to  solve  the  scattering  problem  as  related  to  electromagnetic  radiation.  In 
principle,  this  problem  requires  solving  the  Maxwell  equations,  or  more  precisely, 
the  wave  equation  describing  the  electromagnetic  field,  subject  to  the  boundary  condi¬ 
tions  imposed  by  the  geometry  of  the  scatterer.  The  mathematical  difficulties  are 
extremely  formidable,  and  except  for  extremely  simple  geometries,  there  is  little  hope 
for  obtaining  a  solution. 

Relatively  speaking,  the  case  of  a  plane  electromagnetic  wav  .  incident  upon  a 
homogeneous  spherical  scatterer  represents  a  simple  geometry.  The  German 
physicist  Gustav  Mie  obtained  a  solution  for  this  problem  in  1908.  The  so’ j.tion  is  in 
the  form  of  rather  complex  mathematical  functions,  and  it  is  inconvenient  for  direct 
application  and  ordinary  hand  computations.  The  application  of  the  solution  remained 
relatively  static  until  the  advent  of  modern  high-speed  computers  brought  adequate 
computation  within  the  realm  of  possibility.  Since  the  Mie  solution  is  general,  it  spans 
the  range  of  optical  phenomena  and  properties,  thus  leading  to  rather  diverse  applica¬ 
tions.  Although  several  tables  of  Mie  functions  have  been  generated,  the  use  of  such 
tables  is  reasonably  involved.  Further,  the  coverage  by  any  given  table  is  often  in¬ 
adequate  for  the  specific  needs  of  a  given  investigator. 
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A  Drief  consideration  of  the  Mie  problem  will  demonstrate  why  this  condition 
exists.  In  a  very  strict  sense  the  designation  Mie  functions  should  be  reserved  for  the 
quantities  Sj  {JH>, m,  0)  which,  in  dimensionless  form,  are  the  amplitudes  of  the  farfield 
electric  vectors.  J*  these  quantities  are  divided  by  k  =  271/Xj  where  X  is  the  free-space 
wavelength  of  the  radiation  thus  making  k  the  free-space  propagation  constant,  the 
resulting  quantities  are  the  amplitudes  of  the  electric  vectors: 


Ai  = 


b1  (jm/,  m,  0) 


k  » 


271 

T 


A  - 
2 


S  (>»*',  m,  0) 


i 


Thus,  so  long  as  care  is  exercised  in  the  usage  of  these  quantities  in  a  physical  sense, 
Sj  and  S£  can  be  interpretea  as  electric  amplitudes. 

One  of  the  difficulties  encountered  in  attempting  to  tabulate  these  Mie  functions 
can  be  appreciated  from  the  three  variable  dependence  of  the  quantities.  Actually 
since  the  index  of  refraction  is  a  complex  quan!  .ty,  this  relation  is  more  appropriately 
taken  as  a  four  variable  dependence.  The  total  number  of  values  to  be  tabulated  is  the 
four-term  product  of  the  number  of  values  assigned  to  each  independent  variable.  If 
any  appreciable  range  of  values  is  assigned  to  the  independent  variables,  the  number  of 
values  to  be  tabulated  soon  becomes  unwieldy. 

There  are  several  other  quantities  that  can  be  understood  when  the  description  or 
designation  Mie  functions  are  made.  Two  of  these  are  the  normalized  extinction  and 
scatteiing  cross-sections,  also  termed  the  efficiency  factors.  These  quantities  are  the 
actual  crobs-sections  divided  by  the  geometrical  cross-3ections.  Again  those  two 
quantities  have  a  three  rarameter  dependence  upon  the  Mie  number  and  the  complex 
index  of  refraction,  so  again  a  tabulation  covering  any  considerable  range  for  one  or 
more  of  the  variables  becomes  quite  voluminous  and  even  modest-range  coverage  is 
unwieldy.  The  availability  of  computing  facilities  has  made  the  direct  approach  of 
programming  the  problem  a  more  feasible  approach,  and  this  has  been  so  suggested  in 
the  literature.  The  investigator  can  then  exercise  complete  flexibility  in  studying  his 
specific  case. 


THE  RESEARCH  PROGRAM 


The  objective  cf  this  research  program  as  described  in  the  work  statement  was 
twofold.  The  primary  objective  of  this  study  was  to  provide  quantitative  data  from 
which  first-order  estimates  of  the  reflective  power  of  chemically  induced  clouds  could 
be  made.  A  secondary  objective  was  to  provide  the  means  of  making  first-order 
estimates  related  to  weapons  effects  in  general. 

As  the  program  developed,  the  greater  portion  of  our  efforts  was  devoted  to  the 
primary  objective.  As  a  result  of  the  program,  a  tv.  lputational  procedure  -  a 
computer  program  in  being  -  is  available  for  making  first-order  estimates  cf  the 
reflective  power  of  clouds.  To  date  limited  quantitative  data  have  actua’ly  been  pro¬ 
duced.  The  quantity  cf  data  produced  would  rapidly  increase  were  a  stimulus  for 
such  activities  to  appear. 


. . . 
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In  principle  the  computer  program  has  also  fulfilled  the  secondary  objective.  I*, 
can  be  used  for  weapons  effects  in  general,  but  to  do  so  imposes  upon  the  requesting 
engineer  the  burden  of  expressing,  or  re-expressing,  his  information  in  terms  com¬ 
patible  with  the  computer  program. 


The  Mie  Solution 


The  complete  mathematical  development  of  the  Mie  solution  is  extremely  com 
plex.  Since  the  purpose  here  is  utilitarian,  this  mathematical  development  is  not 
presented  in  detail. 

The  working  equation  describes  the  dpecific  intensity  of  the  monochromatic, 
scattered,  rrdiant  energy: 


P(9) 


where 


Jr  =  incident,  monochromatic  flux,  cal /sec  •  cm^ 
^  =  total  scattering  cross-section,  cm^ 


P(0)  «  angular  dependent  phase  function. 


This  relation  gives  the  amount  of  energy  scattered  by  a  small  volume  element  of  a 
monodisperse  medium  when  the  volume  element  is  illuminated  by  a  parallel  unpolarized 
flux.  Since  the  scattering  is  done  by  a  small  volume  element,  ^)-^{0)  is  expressed  on  a 
per  unit  volume  basis,  and  since  is  dependent  upon'the  scattering  angle,  0,  the 

specific  intensity  is  normalized  by  the  factor  1/4ti  to  a  per-solid-radian  basis: 


^<0>=-^sc-^r  — ^ 


sec 


cm 


strd 


The  total  scattering  cross-section  is  given  by: 

£sc  =  n<r'  asc  r* 


where 


n(r)  =  number  density  of  spherical  particles  having  a  radius  r,  particles/cm^ 

.  •  z 

0sc  =  the  scattering  cross-section,  cm  . 


V 
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The  scattering  cross-aection  can  be  looked  upon  in  the  sense  of  a  fictitious  parameter 
associated  with  a  particle,  or  in  the  sense  of  an  eff  ctive  value  assigned  to  a  param¬ 
eter.  To  illustrate  what  is  meant,  consider  first  an  ordinary  rubber  ball.  As  a 
sphere  it  has  an  unambiguous  value  for  its  diameter.  Having  an  unambiguous  diameter 
means  that  it  also  has  an  unambiguous  cross-sectional  area.  As  the  size  of  the  sphere 
shrinks,  the  meaning  of  the  diameter  is  not  so  unambiguous  and  the  operational 
approach,  that  is,  definitions  that  apply  strictly  only  under  a  specific  set  of  circum¬ 
stances,  must  be  adopted.  For  instance,  an  estimate  of  the  diameter  of  the  particle 
made  on  the  basis  of  an  experiment  concerning  its  volume  may  or  may  not  correlate 
to  the  diameter  determined  by  an  optical  experiment. 

For  a  heterodisper se  media,  the  total  scattering  cross-section  can  be  expressed 
as  a  summation: 


£sc  =  £  n*ri*  asc  <m.  ri>  X)  cm2 
i 


giving 


Jx(0)=^|  ^n(ri)a3c  (m,  rjL,  X,} 


P  (0)  cal/sec 


4tt 


cm* 


strd 


Next,  for  a  given  spherical  particle,  the  phase  function  at  a^v  scattering  angle, 
0,  Is  defined  as: 


P  (**J  m,  0)  =  2 


S1  m,  0) 


S2(**,  m,  0) 


K  {/m,  m) 


where 

Sj  =  dimensionless  amplitude  of  one  farfield  electric  vector  —  bv  con¬ 
vention  Sj  is  perpendicular  to  the  scattering  plane 

S2  =  dimensionless  amplitude  of  the  second  farfield  electric  vector  — 
by  convention  Sg  is  parallel  to  the  scattering  plane.  By  definition 
the  scattering  plane  is  that  plane  containing  both  the  incident  ray 
and  the  scattered  ray;  the  scattering  angle,  0,  is  measured  iu 
this  plane. 

«/*t  =  27ir/X  =  Mie  number 
r  =  radius  ol  spherical  scatterer 

X  =  wavelength  of  the  incident  and  scattered  radiation. 
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The  Mie  number  is  a  generalized,  dimensionlesti  parameter  by  which  radiant  energy 
effects  may  be  correlated.  Since  it  is  essentially  the  ratio  of  the  particle  diameter  to 
the  wavelength,  it  subordinates  both  the  individual  effects  and  emphasizes  the  combined 
effect,  the  scattering  event: 


K8C  =  — —  =  normalized  scattering  cross-section  or  scattering  efficiency 
TTr 


=  ratio  of  scattering  cross-section  to  the  geometrical  cross-section. 


This  I*  a  dimensionless  parameter  which  relates  the  effectiveness  of  the  scattering  to 
the  more  familiar  concept  of  physical  size.  5"nce  Sj  and  Sg  are  complex  quantities,  the 
expression  for  the  phase  function  may  be  writ>en  as: 


Pt*,  m,  0)  = 


Sj  Ut,  m,  m,  0)  f  (*,  m,  0)‘S2*(/m,  m,  0) 

- 

Ksc  Us  m) 


The  phase  function  may  be  looked  upon  as  a  partition  function  in  the  sense  that  it 
delineates  what  portion  of  the  incident  energy  is  scattered  into  a  differential  solid  angle 
characterized  by  a  particular  scattering  angle. 

The  expressions  for  Sj.  S?,  and  KgC  are  infinite  summations: 

00 

51  l**,  n>,  0)  =  Y  {*„  ™a  +  bn  T„} 

n=  1 

oo 

52  t«,  m,  0)  s  V  *  ‘~*r  {b  71  +  a  T  T 

£  ’  n'n+1)  n  n  n  nj 

n=l 

oo 

KSCU.  <n}=-L  ^  (2n+l)|jan|2  +  |bn|2}  . 

n=  1 


As  described  earlier,  the  efforts  here  were  intended  to  provide  a  methodology  by  which 
working  estimates  could  be  quickly  made.  In  line  with  this  objective,  the  computer 
programs  were  written  with  a  maximum  value  of  n  =  10. 

The  mathematical  complexity  increases  sharply  in  the  expressions  for  an  and 
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[  4"  +  S  ]  J„+'/2  *"»  - 


a_  = 


a  r  A 


A  1  r 

—  +  —  J  t 
m  J  |_  n+  * 


y2  (**)  +  i  (-)  J_n_i4  W 


Jn-‘/i  W  +  1  (_)  J-n+1/i  U) 


b  = 
n 


f  mA  +- 
n 


Jn+1/2  M  "  Jn-l/2  ^ 


mA+i  J  ,x/  W+  i  {-)"  J  i 
n  -*»  J  L  n+/2  -n-* 


■/.  w 


Jn-'/2  W  +  1  (-)  J_n+./2  t") 


In  these  expressions 


A 


n 


-n 

m,?n- 


J  ,/  (iTufll) 
n-/a 

Jn+!4 


Trying  to  give  physical  meaning  to  these  coefficient®  would  be  an  extremely  difficult 
job  at  best,  and  to  try  to  do  so  without  a  detailed  derivation  would  be  essentially  im¬ 
possible.  The  wisest  course  seems  to  be  to  merely  accept  these  as  expansion 
coefficients. 


The  expressions  for  71  and  Tn  are  somewhat  simpler  in  form  and  involve  the 
Legendre  polynomial: 


Tn-M*„-  0 -M2)^ 

where 


Pn  (/i)  =  Legendre  polynomial 

H  =  cos  G  - 1  \i  <  1 

Actually,  this  rather  lengthy  discourse  is  really  a  oui'.e  abbreviated  resume  of  the 
mathematics  involved. 


The  Computational  Program 


The  equations  summarized  in  the  previous  section  were  programmed  for  the 
Control  Data  3400  computer.  The  program  is  referred  to  as  PR0GRAM  CL(£UD. 
PROGRAM  CL<JUD  was  written  in  3400  FORTRAN,  a  superset  of  FORTRAN  IV.  Com¬ 
plex  variables  were  used  throughout.  This  program  calls  three  subroutines  BESSELC, 
LECPD1,  and  PLOTXY1.  The  former  two  (for  computation  of  Bessel  functions  and 
Legendre  polynomials)  are  included.  PLOTXYl  should  be  supplied  by  the  user  for  the 
plots  of  J-x(6)  versus  6.  The  general  flow  chart  for  the  program  is  shown  in  Figure  1. 
The  discussion  in  this  section  is  built  around  this  flow  chart  and  provides  a  general 
understanding  of  the  computational  scheme. 
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The  intent  of  this  program  was  to  provide  two  computational  options  referred  to 
as  a  ^ingle-point  computation  and  a  many-point  computation.  In  the  single-point  case, 
attention  is  confined  to  one  value  for  each  of  the  significant  parameters:  radius  of 
scatterer,  scattering  angle,  and  incident  wavelength.  In  the  many-point  case,  at  least 
one  of  these  variables  spans  a  range  of  values. 


The  Input  Parameters 


The  input  parameters  to  this  program  are  the  radius  of  the  spherical  s'atterer, 
the  complex  index  of  refraction  of  the  scatterer,  and  the  scattering  angle.  With  these 
input  parameters  and  for  individual  values  of  the  wavelength  of  the  incident  radiation, 
PROGRAM  CLOUD  computes  the  first  electric  amplitude,  the  second  electric 
amplitude,  the  scattering  efficiency,  the  phase  function,  and  the  scattered  intensity. 

With  the  radius  of  the  particle  and  the  specified  wavelength,  PROGRAM  CLOUD  com¬ 
putes  the  Mie  number,^,  as  a  generalized  parameter  for  discussing  the  interaction 
of  radiant  energy  with  a  scatterer. 

The  input  to  the  program  is  contained  on  three  input  cards,  and  the  presentation 
of  information  on  these  cards  differs  according  to  the  decision  required  from  the  pro¬ 
gram.  The  decision-making  capability  of  PROGRAM  CLOUD  is  based  upon  a  variable 
designated  as  KEY.  ^he  decision  is  raised  in  the  question  "Is  this  case  o  lly  a  single 
point  calculation?  The  answer  "YES"  is  indicated  by  assigning  to  the  variable  KEY 
the  value  of  unity,  KEY  =1.  For  this  case,  the  computations  are  carried  out  for  a 
single  value  of  the  radius  of  the  scatterer,  a  single  value  of  the  scattering  angle,  and  a 
single  value  of  the  wavelength  cf  the  incident  radiation.  The  answer  "NO"  to  the  decision 
question  "Is  this  case  only  a  single  point  calculation?"  is  indicated  by  assigning  to  the 
variable  KEY  the  value  of  zero,  KEY  =  0.  For  this  c?  ie,  computations  are  made  for 
all  combinations  of  all  radii  of  the  spherical  scatterer,  of  all  scattering  angles,  and  of 
all  wavelengths  of  the  incident  radiation. 


The  YES  Case 


For  the  case  where  the  answer  is  YES,  KEY  =  1,  all  of  the  input  data  are  con¬ 
tained  on  the  first  card,  Card  Number  i,  and  the  arrangement  shown  in  Figure  2  is 
used.  The  variables  RER,  RERP,  ALMBT,  RDD,  and  THETT  are  all  floating  point 
variables  written  in  an  E10.  3  format.  The  variable  KEY  is  an  integer  variable  written 
in  an  15  format.  These  format  s^lectio  is  are,  of  course,  quite  arbitrary,  and  selec¬ 
tions  considered  more  appropriate  by  the  individual  users  may  be  readily  substituted. 


The  NO  Case 


For  the  case  where  the  answer  is  NO,  KEY  =  0,  all  three  data  cards  are  needed 
and  the  arra.  ments  shown  in  Figures  3  and  4  are  used.  The  variables  RER,  RERP, 
ALMBT,  RL  ,  and  THETT  are  all  floating  point  variables  written  in  an  E10.  3  format. 
The  variables  NR.  NT,  and  KEY  are  integer  variables  written  in  an  15  format.  As 
before,  the  choice  cf  format  was  arbitrary,  and  appropriate  substitutions  may  be 
made. 


BAT  TEL  LC  20  FIELD 


Columns  1  through  10  -  the  variable  named  RER.  This  Is  the  real  part  of  the  index  of  refraction 
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Columns  71  through  80  -  the  variable  named  ALMBTP.  If  KTAPE  ■  0,  these  colur 
KTAPE  ±  0,  the  value  of  the  wavelength  is  listed  for  which  plots  are.  to  be  made. 
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The  second  card,  Card  Number  2,  is  used  to  introduce  the  values  of  the 
variable  named  RAD,  the  radius  of  the  spherical  scatterer,  and  the  arrangement  is 
shown  in  Figure  4.  These  are  floating  point  variables  so  either  an  E10.  x  or  an 
F10.  x  format  may  be  used.  The  number  of  values  introduced  should  be  e-^ial  to  the 
value  assigned  to  the  variable  named  NR,  up  to  a  maximum  of  twenty-one  values.  A 
maximum  of  three  cards  may  be  used  to  introduce  these  values,  and  the  values  may  be 
distributed  among  these  three  cards  in  any  convenient  manner. 

The  third  card,  Card  Number  3,  is  used  to  introduce  the  values  of  the  variable 
named  THETA,  the  scattering  angle.  The  arrangement  is  the  same  as  for  the 
previous  case  as  illustrated  by  Figure  4. 


The  Output  Parameters 


In  deciding  the  printed  output  for  this  program,  consideration  was  given,  par¬ 
ticularly  to  the  parameters  that  a  user  would  need  in  making  a  physical  application  of 
the  data.  Obviously  the  input  to  each  computation  would  be  needed.  Although,  in 
principle,  the  user  should  have  these  parameters  readily  available  since  he  specified 
the  input  data,  experience  demonstrated  that  having  these  parameters  printed  out  on  the 
data  sheet  by  the  computer  was  a  great  convenience.  Hence,  for  each  pair  of  a  radius 
and  an  angle  of  scattering,  the  following  are  printed:  Radius  of  Scatterer,  Angle  of 
Scattering,  Complex  Refractive  Index,  Initial  Value  of  Wavelength,  Final  Value  of 
Wavelength,  and  the  Increment  for  the  Wavelength. 

Within  the  program,  specific  computation  is  made  of  the  following  parameters: 
Mie  Number,  First  Electric  Vector,  Second  Electric  Vector,  Scattering  Efficiency, 
Phase  Function,  and  Scattered  Intensity.  After  considerable  deliberation  it  was  de¬ 
cided  that  for  the  pragmatic  user,  the  investigator  who  needs  only  a  working  answer, 
only  two  of  these,  the  Mie  Number  and  the  Scattered  Intensity,  are  actually  needed. 
Accordingly,  the  printout  of  results  computed  by  this  program  was  limited  to  these 
two  parameters. 


Plotting  Option 


Examining  many  data  output  tabulations  and  comparing  these  tabulations  are 
rather  burdensome,  and  often  the  comparisons  must  be  plotted  and  studied  further  to 
bring  a  true  understanding.  In  an  effort  to  make  such  comparisons  easier  and  more 
understandable,  a  plotting  op  ion  was  built  into  the  program.  Figure  5  sho-.vs  the  form 
of  the  output  data  as  they  were  plotted  directly  by  the  computer. 
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FIGURE  5.  SCATTERED  INTENSITY  VERSUS  SCATTERING  ANGLE 
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